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ABSTRACT 

As first shown by Shvartsman (1970), a neutron star accreting close to the Edding- 
ton limit must acquire a positive charge in order for electrons and protons to move 
at the same speed. The resulting electrostatic field may contribute to accelerating 
positrons produced near the star surface in conjunction with the radiative force. We 
reconsider the balance between energy gains and losses, including inverse Compton (IC), 
bremsstrahlung and non-radiative scatterings. It is found that, even accounting for IC 
losses only, the maximum positron energy never exceeds ~ 400 keV. The electrostatic 
field alone may produce energies ~ 50 keV at most. We also show that Coulomb col- 
lisions and annihilation with accreting electrons severely limit the number of positrons 
that escape to infinity. 

Subject headings: Acceleration of particles — Accretion, accretion disks — stars: neutron 
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1. Introduction 



2. The Electric Field 



It was suggested a long time ago (Shvartsman 1970; 
Michel 1972; Anile, & Treves 1972; Maraschi, Reina, 
& Treves 1974) that accretion of a hydrogen plasma 
onto a neutron star must necessarily imply the pres- 
ence of an electrostatic field £ under stationary con- 
ditions. The field £ arises because the radiation drag, 
due to the outflowing flux, acts efficiently only on 
electrons while practically all the mass is in protons. 
The electrostatic field compensates the different ac- 
celeration of the two species and allows electrons and 
protons to move at the same speed. Obviously £ is 
stronger the larger is the deviation from free-fall, a 
situation which is expected when the luminosity ap- 
proaches the Eddington limit. 

The typical values of the charge Q associated with 
the electric field are not large enough to influence 
the space-time around the compact object (Michel 
1972), but its presence may nevertheless have interest- 
ing astrophysical consequences. Shvartsman (1970) 
proposed that such a system could act as an electro- 
static accelerator for positrons, for which the electro- 
static and radiative forces have the same direction. 
Balancing energy gains and losses, due mainly to in- 
verse Compton with soft ambient photons, he esti- 
mated that positrons may be accelerated to energies 
of a few tens of MeV. 

Although the "Shvartsman Accelerator" is defl- 
nitely ingenious, it received little attention in the 
past, possibly because no obvious positron source 
seemed available near neutron stars. However, positrons 
may be generated by the radioactive decay of the neu- 
tron star crust elements subject to the bombardment 
of accreting particles (Reina, Tarenghi, & Treves 
1974; Bildsten, Shapiro, & Wasserman 1992). Pair 
production may be also expected in accretion models 
with shocks (see e.g. Shapiro, & Salpeter 1975) or 
in hydrostatic atmospheres kept at mildly relativistic 
temperatures by comptonization (TuroUa, et al. 1994; 
Zane, Turolla, & Treves 1996). 

Here we reconsider Shvartsman's idea, examining 
the role of cooling processes different from inverse 
Compton (bremsstrahlung, non-radiative scattering) 
and the possibility that accelerated positrons survive 
annihilation with atmospheric electrons, escaping to 
infinity. 



Let us consider spherical accretion onto a neutron 
star (M 1.4A'/0, i?,, = 10^' cm) and for simplicity 
suppose that the accreting material is ionized hydro- 
gen. Neglecting general relativistic corrections, the 
emerging luminosity is related to the accretion rate 

by 

i = ,1) 

When L approaches the Eddington limit Le = 
AnGMrripC/aT, the flow dynamics is strongly affected 
by the radiative force (see e.g Maraschi, Reina, & 
Treves 1974, 1978; Miller 1990; Park & Miller 1991; 
Zampieri, Turolla, & Treves 1993, hereafter ZTT). 
Owing to the reduced effective gravity, the flow de- 
viates from free-fall and, for L ~ Le, a settling be- 
haviour appears close to the star surface, where the 
velocity decreases with radius. Even a small devi- 
ation from free fall, which is always present for a 
non-vanishing luminosity, implies that electrons and 
protons are subject to different accelerations so that 
an electrostatic field £ is built up around the star 
(Shvartsman 1970; Zeldovich, fc Novikov 1971; Michel 
1972). Although, in principle, the same argument 
applies to spherical accretion onto a black hole, the 
efficiency is so low (^ 10~* in absence of magnetic 
fields, see e.g. Nobili, Turolla, & Zampieri 1991) that 
departures from free-fall are really negligible for any 
reasonable value of the accretion rate. The expres- 
sion for £ can be easily derived considering the forces 
acting on the two species (e, p) 



and 



GMrUr, 



£e , 



(2) 



(3) 



where Leo denotes the luminosity in the frame co- 
moving with the accretion flow. Equating the two 

accelerations we get 

S ^ ^ 200 ( ^X^V/cm. (4) 

incer^ \ r J Le 

At large radii, where Leo ^ L, the previous expression 
coincides with the result obtained by Maraschi, Reina, 
& Treves (1974, 1978) analyzing the flow equations for 
the two species 



Lax 



exp 



L ii* 
Le r 



(5) 
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An estimate of the largest potential differenee can 
be easily obtained integrating equation (4) assuming 



L = Lt 



(Ay), 



GMr 



200 MV 



(6) 



which shows that, in principle, extremely relativistic 
positrons may be expected to reach infinity. 

3. Compton Scattering 

Here we consider a scenario in which the number 
density of e"*" is not zero at some radius ri„. Positrons 
are immediately accelerated by the electric field and 
gain energy at a rate 



Es = ve£ , 



(7) 



where v is the particle velocity. At the same time, 
positrons arc accelerated by the radiative force and 
suffer radiative losses through inverse Compton (IC) 
scatterings with the radiation bath; e"*"— p and e"*" — e~ 
bremsstrahlung; non-radiative scatterings with the 
flow electrons. In this section we concentrate on IC 
losses, which is the only cooling process taken into 
account by Shvartsman, and we show that, even in 
this simple case, ultra-relativistic energies are not ex- 
pected to be reached. In the next section we will ver- 
ify that bremsstrahlung losses are negligible, but that 
elastic scatterings may be very efficient in transfcring 
energy to ambient electrons and will ultimately stop 
the positron flow. 

The energy lost by a positron for IC scattering with 
a radiation field of energy density U and radiation 
pressure K is 



Eic = -caTl^P''{U + K), 



(8) 



while the gain due to the radiative force is 

Erad = CaTj^P{l + P^)F, (9) 

where F is the total flux divided by c (Gurevich, & 
Rumyantsev 1965, all the radiation moments are here 
evaluated in the lab-frame). For an isotropic radia- 
tion field, F = 0, K = U/3 and the expression for 
the net power reduces to Ejc — — 4c(TtC/7^/3^/3 (see 
e.g. Rybicki, & Lightman 1979). Accordingly, the 
positron energy changes with radius as 



dE 
dr 



= e£ + aTj'^{l+p'^)F-aTiU + K)'y'^(3. (10) 



Introducing the dimensionlcss quantities x = r/vg, 
u = 'lirrgCU/LE, k = 4iTrgCK/LE, and / = L/Le, 
where is the Schwarzschild radius, equation (10) 
can be written as 

dx 2me °° 



(11) 



Whenever x <C x^ = mp/2TOe the solution of equa- 
tion (11) relaxes immediately to the reduced solution 



x'^{u + k) 
2 {Ico - I) 



•1 + 



x^{u + ky 



(12) 



irrespective of the initial particle energy. If no other 
processes prevent positrons from escaping (see the 
discussion in section 4), the terminal Lorentz factor 
can be obtained integrating numerically equation (11) 
once the radiation field in the accreting medium is 
known. In section 5 numerical results based on the 
accretion model described in ZTT will be presented. 

Equation (12) also allows a close comparison be- 
tween our result and Shvartsman's one. The contri- 
bution of electrostatic acceleration can be easily ob- 
tained dropping Erad in equation (10) which is, for- 
mally, equivalent to put Z = in equation (12). An 
upper limit for the terminal Lorentz factor in this 
case is given by 7^ = (l -t- ^/2) /2 ~ 1.2, which 
corresponds to a kinetic energy 50 keV. This is 
much less than the value of tens of McV proposed by 
Shvartsman. The point in Shvartsman's original anal- 
ysis that should be reconsidered is the calculation of 
the IC losses, which were estimated in a very crude 
way by multiplying the energy exchange in each col- 
lision times an "optical depth" parameter, r+. Even 
assuming that the total number of collisions suffered 
by the positron in traveling a radial distance r, N, is 
indeed related to Shvartsman's equation (2) (and this 
is questionable), N should be w r^, since t+ 1. 
This implies that Shvartsman's Compton cooling is 
underestimated by a factor « r+ ~ 10^. Moreover 
the final particle energy is evaluated assuming that 
no losses are present for t-^ < 1 which results in a 
even larger underestimate of IC cooling. 

4. Other Interactions and Relevant Timesca- 
les 

In the previous section we have derived the energy 
equation for positrons, including only IC cooling. In 
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section 5 we show that, cvcm under this favourable 
assumption, the positrons terminal Lorentz factor is 
<^ 2 (see table 1). In order to establish if indeed 
these mildly rclativistic particles may reach infinity, 
we need a more thorough understanding of the phys- 
ical conditions in the flow close to the star surface. 
In our scenario four different kinds of particles are 
present: "primary" electrons and protons of the ac- 
creting material, with the same number density n and 
temperature T, together with positrons and electrons, 
with number density n+ and temperature T_|_, that 
are created in the atmosphere by photon-photon or 
particle collisions. We will limit our discussion to the 
case in which z = n+/n <C 1, so the number den- 
sity of electrons is n_ = n + n+ ~ n. The first is- 
sue that should be addressed concerns the probabil- 
ity that positrons annihilate with ambient electrons. 
This question arises quite naturally since, at variance 
with Shvartsman's scenario in which highly rclativis- 
tic energies are expected, now the particle energy is 
not high enough to ensure that the cross section for 
pair annihilation becomes negligible. 

In the non-relativistic limit, the annihilation rate 
for positrons and electrons with arbitrary distribution 
functions is (Svensson 1982) 



n+ = n+n_C7rrg . 



(13) 



We introduce a characteristic time tA for annihilation 
defined as 



tA = ^ = 

n+ 



1 



n-cnri 



1.34 X 10-3 



10 



17 



n 



(14) 



Since tA depends on n, it gives the probability of an- 
nihilation on a scale of order the density lengthscale 
in the accretion flow, /i r. On the other hand, a 
positron of energy 7 covers a distance /i in a time 



td 



r+h 



7 



-dr ~ 



7 



h, (15) 



r+h/2 



and the particle survives annihilation if td < tA- A 

numerical comparison between these two timcscales 
will be given in the next section for different values of 
both the emergent luminosity and ri„. Here we only 
note that, in the more favourable situations in which 
particles are immediately accelerated up to 7 ~ 1.7, 
this condition is met if 



n < 



1 



irri 



3.2 X 10 



18 



10^ 
r 



cm 



(16) 



This limit is fairly easy to satisfy for near Eddington 
accretion since n = M/(47rr^mpZ;)~ LB/(47rr^mpC"^)« 
10^''(10^/r)^ cm-3 close to the star surface. 

Other cooling mechanisms that could influence the 
maximum positron energy are radiative and non- 
radiative scattering with ions and electrons in the ac- 
cretion flow. The positron-ion bremsstrahlung energy 
loss is given by 



Eb =i -IQ-'^'^ Z'^ncE ergs 



(17) 



where Z is the charge of the ions. This process acts 
on a characteristic timescale 



tB 



E 
Eb 



3.3 X IQ-^Z- 



10 



17 



which, compared with tA, yields 



tA 
tB 



~ 4 X IQ-^Z 



2-7-2 



(18) 



(19) 



so, for pure hydrogen, bremsstrahlung losses are negli- 
gible for particles surviving annihilation and free-free 
could represent a competitive cooling process only for 
large Z. 

On the other hand, the energy exchange due to col- 
lisions is indeed important and can modify substan- 
tially the picture outlined up to now. Let us introduce 
the energy loss rate ve for a test particle in a plasma 
of protons and electrons as 



Eel = —VeE , 



(20) 



and focus on the non-relativistic regime. The more 
efficient process is radiationless, elastic scattering be- 
tween positrons and electrons; positron proton cquipar- 
tition needs a larger timescale by a factor mp/mg. 
The energy loss rate can be obtained from the ex- 
pression of diffusion coefficients, which, in the limit v 
much larger than the proton velocity, yields (see e.g. 
Rosenbluth, & Sagdeev 1983) 



Ve 



87re*nlnA 



me 



-I- erf 



(21) 



where In A is the Coulomb logarithm and Ug the elec- 
tron velocity. The rate i^e is not necessarily positive 
since a fast particle tends to lose energy to the plasma 
and a slow one to gain energy. Initial differences in the 
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mean velocities of positrons and electrons are elim- 
inated by collisions in a time tcoii ~ I'^Bp^- The 
importance of collisions can be quantified comparing 
tcoll with the dynamical timescale given by expression 
(15). 

5. Numerical Estimates 

In order to investigate the effects of electrostatic 
plus radiative acceleration on particles injected in dif- 
ferent atmospheric layers, we have solved equation 
(11) varying both the accretion luminosity and the 
injection radius ri„ . The run of the flow variables has 
been taken from the one-fluid models computed by 
ZTT for near-Eddington accretion of a "cold" , pure 
scattering medium onto neutron stars. The solutions 
are plotted, together with the reduced solution (equa- 
tion 12), in figures la, b, c for L/Le =0.3, 0.6 and 
0.9, respectively. For any given luminosity, we have 
considered three different values of Vm, which corre- 
spond to a scattering optical depth r ~ 0.1, 1, 5 (see 
table 1), with the exception of the case L/Le = 0.3, 
for which the maximum optical depth is r ~ 1.5 at 
the star surface. This simple choice allows us to mimic 
the acceleration of positrons initially created in differ- 
ent atmospheric layers, and, increasing the luminos- 
ity, to visualize the process in the presence of more 
and more expanded atmospheres. In order to account 
for the anisotropy of the radiation field, we introduced 
the closure /(r) = K/U, adopting for / the same ex- 
pression used by ZTT, /(r) = 2/ [3 (l + r^)] -|- 1/3. 
The final energies in all the cases we considered are 
listed in table 1. Figures 2a, b, c show the compar- 
ison between ta and tA for the same choice of the 
parameters. As can be seen from the figures, the re- 
duced solution is always attained at small radii, but 
at the same time annihilation acts efficiently in the 
inner dense atmospheric layers preventing positrons 
created in optically thick regions from escaping. The 
chance to surviving annihilation increases for particles 
created in photospheric regions and above, although 
positrons starting at r <^ 1 reach infinity with a lower 
Lorentz factor. The larger terminal velocities are ex- 
pected for low values of L when the photosphere is 
closer to the star surface. 

Figures 3a, b, c show td together with tcoU for two 
different values of the gas temperature, since the effi- 
ciency of Coulomb collisions strongly depends on the 
difference between the test particle and the electron 
velocities. We note that only for kT ^ 70 keV is 



Ve negative in an extended region near v^n ; ^o the 
possibility that positrons gain additional energy ow- 
ing to collisions appears remote. On the other hand, 
under the more realistic assumption that kT <J 10 
keV, as expected in typical X-ray sources, we find 
that mildly relativistic positrons are "de facto" pre- 
vented from reaching infinity, all differences between 
the mean energies of the two species being washed out 
by the collective effect of soft collisions. 

6. Discussion 

In this paper we have reconsidered acceleration 
of positrons near a neutron star accreting close to 
the Eddington limit. The radiative force and the 
electrostatic thrust due to the non-vanishing electric 
field around the neutron star are opposed by inverse 
Compton losses and elastic scattering with electrons 
of the accreting material. Ultimately, these are the 
main agents responsible for both the acceleration and 
deceleration of positrons. 

The maximiim energy which can be achieved by 
positrons owing to the electrostatic acceleration alone 
is limited to ^ 50 keV which differs substantially from 
Shvartsman's estimate of few tens of MeV. The main 
contribution to the particle's acceleration is therefore 
provided by the radiative force, which acts equally 
on electrons and positrons. We have also shown that 
the largest terminal Lorentz factors are still ^ 2, even 
if this additional gain is accounted for. A first con- 
sequence is that the maximum energies attained by 
positrons are in a range where pair annihilation is 
still important, at least for particles created at r > 1. 
Moreover, elastic scatterings with cold electrons in 
the accretion flow oppose the particle acceleration and 
provide an additional source of trapping. This sug- 
gests that only a fraction of positrons will actually 
leave the star, while the rest should annihilate and 
eventually contribute to a narrow 511 keV line. 

Furthermore, no multiplication through electro- 
magnetic showers either on the accreting material 
or on the photons can be expected, so the positron 
production depends mainly on the injection process. 
Moreover, we note that the same conclusions dis- 
cussed above can be applied to particles injected be- 
low ^ 400i?* (i.e. below Xc), quite independently of 
the initial conditions. If the source is close to the 
Eddington limit, even if positrons are injected with 
energies ^ MeV, like in radioactive decay, because of 
the efficiency of IC losses, their final energy will be 
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limited to ^ 400 1«;V and the possibility of outstream 
or annihilation will not change. 
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Table 1 

The radius, the proton number density and terminal kinetic energy for different values of the 
accretion luminosity and of the scattering optical depth at the injection point 



L/Le 


T 




n{r,,n) 








10" cm 


10-'-' cur-' 


kcV 


0.3 


1.5 


1 


19 


358 


0.3 


1. 


2 


5 


358 


0.3 


0.1 


141 


0.006 


296 


0.6 


5. 


1.5 


61 


293 


0.6 


1. 


13 


0.5 


293 


0.6 


0.1 


855 


0.002 


151 


0.9 


5. 


11 


5 


198 


0.9 


1. 


123 


0.06 


198 


0.9 


0.1 


5000 


0.0002 


50 



7 



REFERENCES 

Anile, A.M., & Treves, A. 1972, Ap&SS, 19, 411 

Bildstcn, L., Salpeter, E.E., & Wasserman, I. 1992, 

ApJ, 384, 143 

Gurevich, L.E, & Rumyantsev, A. A. 1965, Sov. Physics- 
JETP, 20, 1233 

Maraschi, L., Reina, C, & Treves, A. 1974, A&A, 35, 

389 

Maraschi, L., Reina, C, & Treves, A. 1978, A&A, 66, 
99 

Michel, F.C. 1972, Ap&SS, 15, 153 

Miller, G.S. 1990, ApJ, 356, 572 

NobiH, L., TuroUa, R., & Zampieri, L. 1991, ApJ, 383, 
250 

Park, M.-G., & Miller, G.S., 1991, ApJ, 371, 708 

Reina, C., Treves, A., & Tarenghi, M. 1974, A&A, 32, 
317 

Roscnbluth, M.N., & Sagdccv, R.Z., 1983, Handbook 
of Plasma Physics, vol. 1, (Amsterdam: North- 
Holland) 

Rybicki, G.B., & Lightman, A.P., 1979, Radiative 
Processes in Astrophysics, (New York: Wiley) 

Shapiro, S.L., & Salpeter, E.E. 1975, ApJ, 198, 671 

Shvartsman , V.F. 1970, Astrofizika, 6, 309 

Svensson, R. 1982, ApJ, 258, 335 

TuroUa, R., Zampieri, L., Colpi, M., & Treves, A. 
1994, ApJ, 426, L35 

Zampieri, L., Turolla, R., & Treves, A. 1993, ApJ, 

419, 311 (ZTT) 

Zane, S., Turolla, R., & Treves, A. 1996, ApJ, sub- 
mitted 

Zeldovich, Ya.B., & Novikov, I.D. 1971, Relativistic 
Astrophysics, vol.1, (Chicago: The University of 
Chicago Press) 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v4.0. 



8 



Fig. 1. — a) The positron energy vs. radius for 
L/Le = 0.3 and scattering depth at injection 1.5, 
1, and 0.1. The reduced solution is also shown 
(dashed line); b) same as in figure la for L/Le = 0.6 
and r = 5, 1 and 0.1; c) same as in figure lb for 
L/Le = 0.9. 

Fig. 2. The annihilation time Ia (dashed line) and 
the dynamical time td (solid lines) for the same values 
of the luminosity and injection depth as in figure 1. 

Fig. 3. The collisional time t,,oii for kT = 10, 70 kcV 
(dashed lines) and the dynamical time td (solid line) 
for L/Le = 0.3 and r = 1; b) same as in figure 3a for 
L/Le = 0.6; c) same as in figure 3a for L/Le = 0.9. 
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Figure 1a 



Figure 1b 



Figure 1c 




Figure 2a 



Figure 2b 



Figure 2c 
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